Spectroscopic Notation 


Most of the information we have about the universe has come to 
us via emission lines. However, before we can begin to understand 
how emission lines are formed, and what they tell us, we must first 
go over some basic notation and review some atomic physics. We 
will begin with spectroscopic notation. 





Electrons are described with 4 quantum numbers: the principle 
quantum number n, the angular momentum quantum number £, 
the magnetic quantum number, m, and the electron spin, s. No 
two electrons can have the same quantum numbers. 


We can describe the locations of electrons in an atom as being in 
orbitals. Electrons which are in s shell orbitals have zero orbital 
angular momentum. Electrons in p shell orbitals have orbital an- 
gular momentum of 1, electrons in d shells have orbital angular 
momentum of 2, electrons in the f shell have angular momentum 
of 3, etc. The electron’s spin also contributes angular momentum, 
with the values of +1/2. Since angular momentum is a vector 
quantity, the total angular momentum of an atom is the vector 
sum of the angular momenta of all its electrons. 


Only two electrons can reside in an s orbital at once; six electrons 
can be in a p orbital, 10 electrons be in ad orbital, and 14 electrons 
can be in an f orbital. Note that these numbers are not arbitrary. 
In the s orbital, electrons have no orbital angular momentum, 
so only two can fit in (spin-up and spin-down). The unit angular 
momentum of the p orbital allows 6 electrons in, since the angular 
momentum vector can be pointed in any one of three directions (x, 
y, and z). For the d orbital, there are five independent orientations 
for the angular momentum vector (1.e., 2x, 2y, 2z, x+y, and x+z); 
for the f orbital, 7. 





There is a relation between the quantum numbers of an atom. In 
the ground n = 1 state, only s orbitals are allowed; in the n = 2 
state, both s and p orbitals exist; in the n = 3 state, you can s, 
p, and d orbitals. In order, the first few shells of an atom are 


Js 2s 2p 35", 8p s (12.01) 


Rather than specify the quantum number of every electron, it is 
more useful to give a single quantum number for the entire atom. 
Thus, spectroscopic notation uses S, P, D, and F to denote the 
total orbital angular momentum of all the electrons, t.e., 


Lay 4 (12.02) 


When L = 0, the atom is in an S' state; when L = 1, it is in 
a P state, and so on. Note that since angular momentum has 
a direction associated with it, the individual Z values are added 
vectorially. 


Similarly, spectroscopic notation just gives one number for the 
sum of all the electron spins, t.e., 


a (12.03) 


Since electrons can only have spin-up or spin-down, a 2-electron 
system can only have S = 0 or S = 1. 


For light elements, “LS coupling” is a good rule; this simply means 
that the total angular momentum of an atom is the (vector) sum 
of all the orbital and spin angular momenta 


J=L+S (12.04) 


The state of an atom can thus be described using L, J, and S, 
instead of m and s for each electron. In spectroscopic notation, 
we write the state of an atom as 2°+!L,, where 2S + 1 is called 
the “multiplicity”, and L is called the “term”. For instance, an 
atom of neutral carbon in the ground state can be described as 
187,287; 2p? 3 P). 


Since L, S, and J are all quantized vectors, they can have more 
than one orientation in space. In addition, the z-components of 
these vectors are also quantized. This has several implications. 
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Second, for any value of S, 
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| | instance, two electrons with 





the same spin can be in 3 
different states (with 3 differ- 
ent orientations), while two 
electrons with opposite spins 
j have only one possible state. 





Finally, each value of J represents (2J + 1) separate states. Nor- 
mally, since there is no preferred direction for atoms, the energies 
of each of the (2J +1) states are equivalent. However, this degen- 
eracy can be split by giving the atom a preferred direction (say, 
with an electric or magnetic field). 


Note that when a shell is filled, the electrons occupy all possible 
orientations, so that the net L and © of the shell is zero. Thus, 
only unfilled shells contribute to the L and S of the atom. Simi- 
larly, from symmetry, a shell with all the orbitals filled except one 
looks and acts the same as a shell with only 1 (anti-) electron, and 
shells that are missing two electrons are equivalent to shells that 
contain two (anti-) electrons. 


As an example, consider a 3P term. Since it is a P term, L = 1, 
and, from the multiplicity, S = 1. It therefore contains (2L + 
1)(25 + 1) = 9 separate states. Since L = 1 and S = 1 can add 
vectorially to J = 2,1, or 0, these states are ?P2, which contains 
2J +1, or 5 of the states, ?P,, which holds three of the states, and 
3 Py, which is a single state. The value w = (2J +1) is called the 
statistical weight of a level, since it gives the number of degenerate 
states at a single energy level. 


Transition Rates 


The probability of an electron in an atom making a transition from 
state n to state n’ is given by the Einstein probability coefficients. 
For n > n’, these are 


Ann’: the spontaneous emission coefficient 
Bnn: the induced emission coefficient 
Bnm: the absorption coefficient 


In other words, the Einstein A value is the probability of an elec- 
tron decaying from state n to state n’ all by itself. The Einstein 
Bnn’ value deals with stimulated emission: a passing photon with 
an energy equivalent to the jump from n to n’ can prompt (i.e., 
stimulate) an electron decay. (This is the way lasers and masers 
work.) Finally, the B,-, values deal with the probability of ab- 
sorption. All three coefficients are fundamental constants asso- 
ciated with atomic levels; Ann’ has the units of inverse seconds. 
These coefficients are related: it can be shown (with relatively 
little effort), that 
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where w is the statistical weight of the level. High A values 
(such as 108 sec!) are permitted transitions; on average, if A = 
108 sec~+, it will take an electron 1078 seconds to decay. Also, in 
general, the larger the energy difference of a permitted transition, 
the larger the A value, with A « v’. 





Note that not all transitions in an atom are equally probable. 
From quantum mechanics, the most probable transitions must 
have An 4 0, A€+1, and Aj+1 or 0 (unless j is zero to start with, 
in which case j must change by +1). There are several reasons 
for these rules, including conservation of angular momentum (the 
photon carries away 1 unit of j), and the ability of the atom to 
radiate via dipole radiation. If these rules are not satisfied, the 
transition is classically “forbidden”. In practice, that means that 
the A value for the transition is very low, typically AŚ 1 sec. 








Aside: most astronomers who analyze the absorption lines in stel- 
lar atmospheres do not work with the Einstein coefficients per se: 
instead, they perfer to describe the strengths of atomic transitions 
via the line’s oscillator strength, f. The two are related by 
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where e is the charge of the electron and €o is the permittivity of 
free space (€9 = 1 in cgs units). 


The Photo-ionization Cross Section 


To understand the physics of an emission region, one needs to 
know two functions. The first is the hydrogen atom cross section 
for photo-ionization from the ground state. It takes 13.6 eV of 
energy to ionize hydrogen; this corresponds to a wavelength of 
912 A, and a frequency of 3.29 x 10!° Hz. At this energy, the 
atom’s cross section to a photon is Ag = 6.30 x 10718 cm?. At 
higher energies, the cross section decreases, roughly as 1/v°, so 
that for hydrogenic atoms, 


i= (3) a (12.08) 


Wavelength (Å) 
100 300 1000 3000 10000 


Relative Flux (ergs/cm?/s/Hz) 


16.5 16.0 15.5 15.0 
Log Frequency (Hz) 


14.5 14.0 

The implication is that extreme-UV light with wavelengths just 
shortward of 912 A are absorbed very rapidly in the neutral hy- 
drogen in the interstellar medium. Consequently, at these wave- 
lengths, only a few objects in the sky are visible. 


Because of v~? dependence of absorption, at shorter wavelengths 
(say, in the soft X-ray band), absorption by interstellar gas is much 
less important (although it still cannot be ignored); by the time 
you reach the hard X-ray and y-ray bands, the ISM is virtually 
transparent. Main-sequence stars have temperatures such that 
only the very tail of their blackbody curves extends shortward 
of 13.6 eV. Consequently, the “average” cross section for photo- 
ionization in these regions is quite large. 





The Recombination Coefficient 


The other parameter that is needed is the recombination coeffi- 
cient for hydrogen, a. The rate at which free electrons will re- 
combine into (any atomic level of) hydrogen depends on how fast 
the electrons are moving. 


Fast moving will sweep 











out a larger volume of = a 7 
space than slow moving 4 (| lo) 
electrons, SO A X VOnL.- 





On the other hand, the cross sectional area for recombination, 
Ont must be inversely dependent on velocity, due to electrostatic 
focusing of the electrons by the protons. The calculation of this 
effect is straightforward: 














An electron that starts out moving towards a proton with impact 
parameter, b, and velocity, v, feels a force perpendicular to its 
motion, 


dv eo. eb eb 
BN ge gO (12.09) 


Integrating this once over dt = dx/v gives the velocity perpendic- 
ular to the direction of original motion, 
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Integrating the equation a second time gives the difference be- 
tween the initial impact parameter, and the impact parameter at 
closest approach 
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Thus, the change in the cross sectional area is 

Ao = Tb’ — n (b — Ab)* = 2rbAb — rAd? (12:12) 


which, in the limit of large impact parameters and small deflec- 
tions, means that 
Ao = 2nbAb x v~? (12.13) 


Thus, due to electrostatic focusing, the cross section for recom- 
bination is inversely proportional to the square of the electron 
velocity. 
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The two effects combined imply that a x v/v? « v~+. In practice, 
not all electrons have the same velocity: the distribution of veloc- 
ities is Maxwellian. (It takes electrons in the interstellar medium 
only a few hours to become totally thermalized.) Since the charac- 
teristic temperature of the electrons, Te x v7, this means that the 
recombination coefficient, a « TY E Thus, a is a weakly depen- 
dent on temperature. Numerically, the recombination coefficient 
for all levels of the hydrogen atom is 


aa = 6.82 x 1071? cm? s7! for T, = 5,000 K 
aa = 4.18 x 107+? cm? s7} for T, = 10,000 K 
aa = 2.51 x 107+ cm? s71} for Te = 20,000 K 
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Parameters of a Stromgren Sphere 


Let’s assume that we have a 40,000 K B-main sequence star sur- 
rounded by an interstellar medium composed entirely of hydrogen 
(density N(H) = 1 particles cm~°). 


How long will a neutral hydrogen atom stay neutral? The photon 
luminosity of a B main sequence star is 
L 
/ — dv = 5 x 10*° photons s7! (12.14) 
vo. Pw 
where Vo is the ionization frequency for hydrogen. Let’s now as- 
sume all the ionizing radiation comes from the star. The photo- 


ionization rate at a distance r from the star is 


P = N(B?) f = d (12.15) 
= —~—a, dy 
vp aar2hy ~ 

Note the terms in the equation. The photoionization rate is pro- 
portional to the number of neutral atoms that are around to be 
photoionized, the number of photons available for photoionization, 
and the cross section for ionization. For a precise calculation, we 
should integrate this equation, but for now, we can just adopt an 
“average” value of a,. For a 40,000 K star, only the high-end 
tail of the blackbody function extends shortward of 912 A, hence 
almost all ionizations occur with a, © 5 x 10718 cm?. So 

L aoe ame 


P = N(H° — idv = N(H°® a, f at 
( if Anr2hy Y eve vo hw 4nr? 
(12.16) 


If we now pick a position, say, 5 pc away from the central star, 
then we can just plug in the numbers to get P ~ 1078 ioniza- 
tions cm~? s71. In other words, in a cubic centimeter of space, 
the time between ionizations is 10° seconds: a neutral hydrogen 


atom can stay neutral for over a year! 
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What is the distribution of states for hydrogen? The time it takes 
an atom to decay from state nL to state n’L’ is just one over the 
Einstein A value. The time it takes at atom to decay from state 
nL to any state is therefore 


1 
L ` `“ Anl nib 


ni<n L'=L+1 


(12.17) 





(Note that by specifying that L’ = L +1, we are neglecting the 
forbidden transitions.). Typical A values for permitted transitions 
are 10*S Anz nr Ñ 108 sec7!, so 10742 Trp 1078 seconds. 


There is one exception to this. A hydrogen atom in the 2 7S 
state, has no permitted way to decay down to 17S. Atoms in this 
state must either be a) collided out of the state, or b) decay via a 
low probability (A = 8.23 sec~') two-photon emission, where the 
atom makes a temporary state for itself between 17S and 27S. If 
this occurs, the atom will be trapped in 2 7S for T = 0.12 seconds. 


In any event, the time for a hydrogen atom to decay is less than 
a second, whereas it must wait a year before a neutral atom gets 
photoionized. Moreover, in the low-density region of the inter- 
stellar medium, the timescale for collisions between electrons and 
atoms (in order to excite an electron to a higher energy level) is 
much longer than the timescale to decay. As a result, essentially 
all the atoms will be in their ground state. The ISM is not in 
thermodynamic equilibrium! 








Finally, note that the energy difference between the ground state 
and first excited state of hydrogen (10.2 eV, or A = 1216 A) 
corresponds to 

he 1 


a VY 5k = T ~ 80,000 K (12.18) 
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This is much hotter than the electrons in an H II region (which, 
as we will see, typically have temperatures of T ~ 10,000 K), so 
collisional excitation of hydrogen from n = 1 to n = 2 is unim- 
portant. Again, this means that all the neutral hydrogen will be 
in the ground state. 


What fraction of hydrogen atoms will be neutral? Let € be the 
fraction of neutral hydrogen, i.e., € = N(H°)/N(H), and again, 
let’s consider a spot 5 pc from a 40,000 K main sequence star. As 
we have seen, at that location 


P=N(H°) | 


Vo 


CO 


„ dv = N(H?) 1078 sec7! 12.1 
Triho YY ue (12.19) 


For an emission region to be in equilibrium, the number of ioniza- 
tions (in each cubic centimeter) must equal the number recombi- 
nations. The latter number is given by 


R=N.Npaa (12.20) 


In other words, the recombination rate is proportional to the 
density of free electrons, the density of hydrogen atoms need- 
ing an electron, and the recombination coefficient, a4 ~ 4 x 
10718 em? s~!. Ionization balance implies 





E 
0 
N(H J TR” dv = NeNpaa (12.21) 
In a pure hydrogen nebula, Ne = Np, and N, = (1 — €)N(#). 
Thus 
€N(H)(107*°) = (1-8 N(HY (4x 107°) (12.22) 


or 


2.5 x 107-*€ = (1 — £P N (H) (12.23) 


For a typical ISM value of N(H) ~ 1, € = 4 x 107°. Virtually all 
the hydrogen is ionized! 
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How thick is the transition region between the ionized and neutral 
material? The mean free path of an ionizing photon is given by 

c= : (12.24) 

= N(A°)a, l 

(The larger the cross section, or higher the density of neutral 
atoms, the shorter the path.) By definition, in the transition 
region, half of the hydrogen is neutral. Thus the mean free path 
of an ionizing photon in this region is 


j 1 
trans “" 0.5N(H)ay 





(12.25) 


For ay ~ 5 x 10718 cm? and N(H) ~ 1, lrans S 4X 1017 cm, 
or 0.1 pc. Thus, when the material starts to become neutral, the 
ionizing photons are rapidly eaten up. Considering that the entire 
H II region might be ~ 10 pc across, the transition region is very 
thin. 


(Note: if the source of ionizing radiation is very hard (for instance, 
as in the power law continuum of an AGN), then the effective value 
of a, can be significantly smaller. This can translate into a large 
transition region.) 





What is the size of an ionized region? To calculate the size of an 
emission region, let’s start with the full equation for ionization 
balance. The ionization rate at every location in the nebula must 
be equal to the recombination rate at that location. 





© And, 
nH) f ~ a, dv = Ne Npa A (12.26) 


(0) 


Here, J, is the mean intensity of radiation, in units of energy per 
unit area, per unit time, per unit frequency, per unit solid angle. 
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(The 47 term is the implicit integration over all directions.) Note 
that J, has two components. Ionizing radiation can come from 
the central star, but it can also come from a diffuse component. If, 
somewhere in the nebula, an electron recombines from the contin- 
uum directly into the ground state, the photon emitted will have 
enough energy to ionize another atom. 


Keeping track of the diffuse ionizing field is fairly complicated. 
However, for many nebulae, there is an excellent approximation. 
Since the cross section for recombination is inversely proportional 
to velocity, most recombinations will be from slow moving elec- 
trons. Those recombinations that go directly into the ground state 
will therefore produce a photon with v ~% vo; such a photon has a 
very high cross section for re-absorption. Therefore, let us make 
the on-the spot assumption that all photons produced by 
ground state recombinations are absorbed locally. Any 
recombination to the ground state produces a photon that is im- 
mediately re-absorbed by a neighboring atom: it’s as if no ground 
state recombination ever occurred. 





(In low density nebulae with N(H) ~ 1, this approximation is 

not very good, since, as we have seen, the mean free path for an 
ionizing photon is 

1 

f= ——__ 12.24 

N(H?)a, ( ) 

If the neutral fraction is small, € ~ 1074, then @ can be large. 

However, as the density of particles increases, the mean free path 


decreases, and the approximation becomes better and better.) 





If we neglect the diffuse component of the ionizing radiation field, 
then the only ionizing photons are those that come from the 
central star. If each square centimeter of the star produces a 
monochromatic flux of rF, (ergs cm~? s-t Hz~*), then any point 
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in the nebula sees this field, diluted by a geometrical term, and 
by the attenuation that has already occurred. In other words, 


OO 4 ; OO F, 2 
nH) | ce dv= N(R°) | 2 (5) aye ™™ dv 


hv hv \ r? 

(12.27) 
where R is the radius of the star, r is the distance to the star, and 
T, the optical depth, reflects the fact that some photons from the 
star have already been absorbed before reaching r. Specifically, 
the equation for the optical depth is 








O VO 





Ae f NU Crd (12.28) 
0 
Since the luminosity of the star is 
L, = 4nR?(nF,) (12.29) 


the equation of ionization balance simplifies to 





N(H° OP Ts, 
l f —aye ™ dv = NeNpaB (12.30) 
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where now we are using ag, which is the recombination rate to 
all levels except the ground state of hydrogen. 


Now, assume that the density of the emission region is roughly 
constant. We can re-write (12.30) as 


oO is 
N(H°) / Ape dv = 4nr*N.Npap (12.31) 


O 


and integrate both sides over r, with R representing the extent of 
the nebula. 


R © T, R 
J NH) f T, e ™ dvdr = f Anr’ NeNpag dr (12.32) 
0 Vo 0 
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We can replace dr with N since 
j d 
wr) =f Nadr = = NCH) a, (12.33) 
0 r 


At the edge of the nebula r = oo (by definition), so 


R oo oo oo 
L dt L 
N(H® aye 7’ ———__ = K —Ty 
| ( if hv” : N(A° ja, a hv vf © ‘ 


(12.34) 

or just 

E 

—dv = Q(H® 12.35 
[a= oH”) (12.35) 
In other words, the left hand side of (12.32) is just the number 
of ionizing photons coming from the star. (This should not be 
too surprising, since, one way or the other, all the photons must 
be eaten up.) In the case of a constant density nebula, the right 
hand side of (12.32) is straightforward to integrate, and, since the 


ionization fraction inside the nebula is extremely high, Np ~ Ne ~ 
N(H). So 


R 
4 
Q(H°) =| N-NpapdV = z1R N(H) ag (12.36) 


Thus, we have a very simple expression for the size of the nebula. 
Our 40,000 B-main sequence star will generate an H II region 
~ 50 pc in radius. 
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Modifications Due to Helium 





Photons with energies greater than 24.6 eV (504 A) can ionize 
helium and, due to the v~? dependence on cross section, such 
photons are ~ 8 times more likely to ionize helium than hydro- 
gen. At first glance, this would seem to mean that helium can 
significantly reduce the size of a Stromgren sphere. However, the 
energy difference between the n = 1 state of helium, and the n = 2 
state is about 19.9 eV. (The exact difference depends on whether 
we are talking about singlet states, or triplet states, and whether 
the n = 2 term is a P or an S.) Thus, almost every recombination 
of a helium atom will produce another photon capable of ionizing 
hydrogen. Consequently, the effect of helium on the size of the 
Stromgren sphere is minor. 














(Note: some interesting physics is associated with the He I atom. 
Since there are two electrons, He I has singlet levels (opposite spin 
electrons) and triplet levels (same spin electrons). Like hydrogen, 
singlet decays from 2 tP to 11S are no problem, but 21S to 115S 
is forbidden; these atoms can only decay by creating a temporary 
state somewhere between 2 +S and 11S and emitting 2 photons 
(A = 51 sec™t). When this occurs, 56% of the time, one of the 
photons will have an energy greater than 13.6 eV. Atoms in the 
triplet state, however, have a harder problem: once they get to 
2 3S, the only way for them to decay is to have one electron 
change its spin. The A value for this is A = 1.27 x 1074 sec™!. If 
the spin-change does occur, then the resulting photon can ionize 
hydrogen, but since the timescale for this is long (7 ~% 2.2 hours), 
it’s possible that the atom will get collided out of 2 35S before it 
can decay. In fact, some of these collisions will dump the electron 
directly in 1 1S. If this occurs, an ionizing photon will be lost. 
Depending on the density, He I will eat up between 4% and 34% 
of hydrogen’s ionizing photons. ) 
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The accurate computation of ionization balance for a hydrogen- 
helium nebula can be moderately complicated, but the problem 
can be approximated by noting that a,(He®) is always signifi- 
cantly larger than a,(H). Thus, if you make the assumption that 
all photons capable of ionizing helium do ionize helium, then 


Co i 4 
[ee = QHe) = Fa (Het Nean(He Ries (1237) 


1 
where vı is the ionization frequency of helium, and Rye+ is the 
radius of the He I Stromgren sphere. 


If the central star is hot enough to produce photons with energies 
greater than 54.4 eV (228 A), then helium can become doubly 
ionized. Because of the extremely large differences between the 
energy levels of He II, a recombination to He II will produce, on 
average, 1.25 photons capable of ionizing hydrogen. So, while a 
small percentage of ionizing photons may be lost due to collision 
processes in He I, another small percentage can be recovered with 
He II. 
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Photo-ionization of Heavy Elements 


An insignificant amount of the central star’s photons go into ion- 
izing the heavy elements. The photo-ionization cross section of 
these elements is complicated, since 1) the experimental values 
for the cross sections are not well determined, 2) the cross sec- 
tions are not smooth functions, due to resonances, 3) other ef- 
fects can place the electrons of metals into excited states, from 
which the ionization can occur. Metal recombination rates are 
similarly uncertain, due to effects such as dielectronic recombina- 
tion. (Imagine a recombination where the photon produced via 
the recombination goes directly into exciting one of the ground- 
state electrons. If this happens, the atom can have two excited 
electrons at once, with a total excitation energy greater than the 
ionization energy. If this occurs, one electron may give its energy 
to the other, causing the atom to re-ionize itself!) 





Below are some of the energies (in eV) needed to ionize metals: 


Element II III IV V 


Hydrogen 13.6 

Helium 24.6 54.4 

Carbon 11.3 244 47.9 64.5 
Nitrogen 145 29.6 47.4 77.7 
Oxygen 13.6 35.1 54.9 77.4 
Neon 21.6 41.0 63.5 97.1 





Note that to strip one electron off oxygen requires 13.6 eV — 
exactly the same as hydrogen. Thus, when hydrogen is ionized, 
oxygen is ionized. 
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Nebular Heating and Thermal Equilibrium 





Consider the energy contained in the free electrons of a nebula. 
Every ionization creates a photo-electron with energy 


1 
5 Me 


= h(v — vo) (13.01) 





while every recombination removes FMeVe. The difference be- 
tween the mean energy of a photoionized electron, and that of a 
recombination electron must be lost by radiation if the nebula is 


to be in thermal equilibrium. 





There is only one source of heating for a nebula: the photo- 
ionization of electrons. It is not difficult to estimate the initial 
temperature (T;) of these electrons. For a pure hydrogen neb- 
ula ionized by a blackbody the average energy of a photoionized 


electron is 20 y 
/ — h(v — vy) dv 
w tw 
kh = se mMm (13.02) 
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One can find an analytic expression for this, if one assumes that 
the ionizing spectrum is a blackbody, and all the ionizing photons 
are on the Wien part of the blackbody curve, i.e., kT, < hvo. In 
this case, 





2 
2 


3 3 


a) ee | o (13.03) 


If you make this substitution, and crank through the math, then 
you find that 


N 


SAT. x kh > Tx zT» (13.04) 
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Two things to note. First, T; depends only the shape of the ion- 
izing radiation, not on its strength. (In other words, the thermal 
energy of an ionized electron depends only on the energy of the 
photon that hit it, not on how many photons there are.) Thus, 
the temperature of a nebula being ionized by 100 large 30,000 K 
stars will be less than that ionized by one 100,000 white dwarf. 
(The size of the nebula will be larger, but the temperature will be 
smaller). 


Second, because the cross section to photoionization decreases as 
v—?, high energy photons will penetrate deeper into a nebula than 
low energy photons. Thus, 7; may increase as the distance from 
the star increases. 





The exact heating of a nebula is not difficult to compute. If G(H) 
is the energy gained by the nebula through the photoionization of 
hydrogen, then 


G(H) = N(H°) A a h(v — vo)ay dv (13.05) 


Vo 





which, through photoionization equilibrium, is 


Ahy — V9 )a, dv 


G(H) = NeNpaa }, 


oo (13.06) 
~ hv” 
Since the last term is simply injected energy, 

G(H) = N.Npaa (Sen) ~x NeNpaa kT, (13.07) 


A similar equation holds for helium, so the total heating is 

G = G(H) + G(He°) + G(He’) (13.08) 
The amount of metals in the nebula is so little, that heating from 
those elements can easily be neglected. 
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Sources of Cooling 


There are three different mechanisms a nebula has for cooling 
itself. The first is through thermal bremsstrahlung of the free 
electrons; the full expression for this is 


25 are® Z2 ( Irk 


1/2 
= 1/2 
ESP alae ) Tel of ¢NeN+ 


Me 
= 1.42 x 10-27 ZT? gff NeN, ergs cm? s! (13.09) 


where gy is the gaunt factor. Since virtually all of these photons 
are longward of vo (most are in the radio), all this free-free lumi- 
nosity leaves the nebula. However, this energy is mostly incon- 
sequential. Free-free emission is the least important of all energy 
loss mechanisms. 





A second source of cooling comes from the recombination of elec- 
trons into hydrogen (and helium). This is computed simply from 
the kinetic-energy weighted recombination coefficient to each level 
of the atom 





Lr(H) = NeNp kT Ba (13.10) 

where i 
= Y ar (13.11) 

n=1 L=0 


A typical value for 6 is 10718 cm? s~!. Note that even through 8 
represents the energy released by a recombined electron, it does 
not have the units of energy. By definition, the electron energy 
has been divided by kT, which is why (13.10) has a kT in the 


expression. 


Recombination cooling is significant, but it does have a very inter- 
esting property. Recall from our discussion about recombination 
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(12.13) that, due to electrostatic focusing, slow moving electrons 
are more likely to recombine than fast moving electrons. Thus, re- 
combinational cooling will continually remove electrons from the 
low-energy tail of the Maxwellian distribution. Consequently, if 
this were the only source of nebular cooling, the temperature of 
the free electrons would actually increase! 


This same effect makes the on-the-spot approximation for thermal 
balance less accurate than it is for ionization balance. (The slow 
moving electrons from the diffuse field are more important.) Still, 
one can use pg, instead of 6,4 in (13.10) and get a reasonable 
approximation for recombinational cooling. 
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Collisional Cooling 


The most important physical process that occurs in an emission 
line region is collisional cooling. Consider an electron in some 
atom in state nL. If nL is not the ground state, that electron 
may decay, at a rate determined by the A value of the transition. 
If nL is the ground state, then, of course, the electron will just 
sit there. However, if the density of particles is high, or if the A 
value for downward transitions are low, a free electron may come 
by and collide the atom out of state nL into state n’L’ before a 
spontaneous transition can occur. If state n’L’ has less energy 
than state nL, then the free electron must carry away the energy. 
However, if state n’L’ has AE more energy than state nL (as it 
must, if the bound electron is already in its ground state), then the 
colliding electron will lose that amount of energy in the collision. 
The nebula will have therefore cooled. 





The point to notice about this effect is where it occurs. The energy 
difference between the ground state of hydrogen and the first ex- 
cited state is AF = 10.2 eV. Therefore, for collisional cooling to be 
important for hydrogen, kT ~ 10.2 eV, or T ~ 120,000 K. Since 
this is much hotter than the temperature of the free electrons, col- 
lisions with neutral hydrogen atoms (which are, of course, rare) 
cannot cause an excitation and therefore cannot cool the nebula. 








For helium, the energy difference between the ground state and 
the n = 2 state is even higher (19.8 ev), so likewise, ground-state 
helium collisions are unimportant. Collisional cooling on metals, 
however, is another story. Ions with electrons in partially-filled 
p-orbitals (such as oxygen, nitrogen, carbon, neon, argon, and 
sulfur) all have energy levels that are ~ kT above their ground 
state. Thus, they are an extremely efficient cooling mechanism. 
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To understand how efficient these ions are, we must consider the 
rate at which these collisions occur. The computation of the rates 
has four terms: 


a term which includes the density of electrons (Ne), and the den- 
sity of ions (N;). 


a term which includes the volume being swept out by a free elec- 
tron each second (o « v) and the effects of electrostatic focusing 
(o x v7). In terms of temperature, this means that the rate of 
collisions must go as 

gT i? (13.12) 


a term which includes the Boltzmann factor. For a free electron 
with velocity v, if AE > m,.v?/2, then clearly, there is not enough 
energy in the collision to excite the electron. In terms of tempera- 
ture, the larger the value of AF, the fewer the number of electrons 
that can cause an excitation. On the other hand, if the collision 
involves a downward transition, there is no energy barrier. Thus 








Gee IP APSO (13.13) 


a term which involves the statistical weights. If level n’L’ is actu- 
ally composed of 10 different states, then the collisional rate out 
of that level will be 10 times smaller that to a similar level with 
only one state. In other words, q x 1/wi. 


a (dimensionless) quantum mechanical cross section, which de- 
pends on the wave functions of the states, q x Q(i, 7). 
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If we put these terms together, the rate of collisions from state 2 
to state 7 is 
Re = NeNiqij (13.14) 


where 











( on yo h2 Qi, j) 
qij = 





kT. m3/? Wi 
ee 
= 8.629 x 1078 G i) cm? s7! (13.15) 
/2 
wk Te 


if state j has lower energy than state i (i.e., AE < 0). If, on the 
other hand, AE > 0, then 


Oli; 
qi; = 8.629 x 107° u EAR ene B (13.16) 
wile 





Interestingly, the quantum mechanical collision strength, Q, is 
symmetrical, i.e., Q(i, 7) = Q(j, i). This causes the two collisional 
rates to be related 


Wi a 
qij = —e Seas (13.17) 


Also, if either S = 0 or L = 0, then there is a simple relation for 
the collision strength between a term with a singlet level and a 
term containing several levels, t.e., 





GIAI) 


ET) E 


Q(SL,S'L') (13.18) 
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For example, let’s look at the energy levels of doubly ionized oxy- 
gen, The ground state of OTT is a “p-squared” configuration, i.e., 
the electrons are in 1s?2s?2p°. In the 2p” shell are 5 different en- 
ergy levels: the ground state term is at L = 1, and, depending 
on the spins of the electrons, the term has three different energy 
states, 3 P), 3P}, and 3P2. A couple of eV above these P states is 
a ! Də state, and about 3 eV above that is a singlet L = 0 state, 
1 Sọ. All these states have n = 2, so transitions between them are 
classically forbidden; in practice, their A values are < 0.1 sec~'. 








Because the D state and the S state are both singlets, a single 
number can be given for the collision strength between °P and 
1D: in this case, O(?P,' D) = 2.17. The collision strength from 
the individual levels are therefore 





QP, D) = =00P} D) 


3 
QË P, D) = 5 2CP, D) 
1 
QË Po,’ D) = 5 2CP,' D) (13.19) 
If the density of particles is low, then it is unlikely that a collisional 
de-excitation will occur before a spontaneous decay occurs. In this 


case, every collisional excitation creates a photon that escapes the 
nebula. The total collisional cooling for a given species is then 


Lo = NO NeNig1,3 hvi,j (13.20) 


J 





where the sum is taken over all excited levels (typically 4), and 
hvı j is the energy difference between the ion’s ground state and 
state 7. 
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Collisional Cooling and Detailed Balance 


If every collisional excitation from the ground state to level 7 were 
followed (eventually) by the emission of a photon (or multiple 
photons) which leave the nebula, then the collisional cooling rate 
would be 





Lo = NeNign; : hiy; (13.21) 


But, in practice, collisional de-excitation also occurs. Consider a 
two-level atom of an ion. Let N; be the number of ions in the 
ground state, and No be the number of ions in the excited state. 
In equilibrium, the number of electrons going out of N; is equal 
to the number of electrons entering level N4, either by collisions 
or by radiative decay. In other words, 


NeNiqi2 = Ne Noq21 + NoAa1 (13.22) 


The relative level population is therefore 


N. Ne 
(=) B ae (13.23) 
Ny Neqai + A21 
The collisional cooling rate for this two-level ion is therefore 
Ni Ne 
Lo = Nə Ao h Yo, = OTTA Ao h V21 (13.24) 


Let’s look at this equation a bit more closely. If we divide both 
the numerator and denominator of (13.24) by Ag, then 


Neq21 


1 
Lo = Ne Nıqı2 h VIS w (13:25) 
A21 TA 


Now note the limits. As the electron density, Ne — 0, collisional 
de-excitation becomes unimportant, and the cooling rate becomes 


Lo — NeNiqi2 hnis (13.26) 
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In this case, the cooling is proportional to the electron density, and 
every collision upward creates a photon which cools the nebula. 
On the other hand, as Ne —> ov, 


Lo —> Ni (=) Ao ih Və] = Ni (=) g OEIED Asih V21 


q21 Wy] 
(13.27) 
This is nothing more than the cooling rate for a gas in thermody- 
namic equilibrium. 


Of course, although a few atoms (such as Be II) can be treated 
with two levels, most require more. In fact, the p?, p°, and pł 
configurations all have 5 low-lying levels, separated with two ma- 
jor divisions. To calculate the cooling from an n-level atom, you 
have to solve a set of linear equations, equating the number of 
electrons entering a level to the number of electrons exiting the 
level. 


» NeNjQji + Ss Nj Aji = D NeNiqij + ` N;Aiz (13.28) 


j#t j>t j#t j<i 


The first sum is the number of electrons entering level į through 
collisions from all the other levels. The second term is the number 
of electrons entering level ¿ through radiative decays from the 
levels above 27. On the other side of the equal sign is the number of 
electrons exiting level 2 via collisions, and the number of electrons 
exiting level ¿ via radiative decays downward. This is the equation 
of detailed balance for a collisionally excited line. 





Note that there is actually one other equation to detailed balance 
S N=N (13.29) 
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This simply says that the sum of the ions in all the levels must 
add up to the total number of ions there are. So this gives you 
n+ 1 equations, and n unknowns. Thus, with some algebra (and 
a little bit of patience) you can solve for the number of ions in 
each level. 


Once the ion densities are known, the collisional cooling follows 
simply from 


Lo = >, N; X Aijhvij (13.30) 


j<i 


The equation of detailed balance also allows us to specify when 
collisions out of a level are more important than radiative transfers 
out of a level. From the left side of (13.28), this occurs when 


`“ NeNjqji > DD Nj Aji (13.31) 


jži j<i 


or, more to the point, 
3 Aij 
i lis 


j#t 


Ne (crit) 





(13.32) 


N- (crit) is called the critical density for the level. If Ne < Ne(crit), 
then collisional de-excitation is not important. For most forbidden 
lines, this occurs, very roughly, at Ne ~ 10° cm7?. 
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Thermal Equilibrium 
A nebula in thermal equilibrium has 
G=Lrp+Lr+Le (13.33) 


where G is the heating due to photoionization (equations 13.06 
and 13.08), Lys is the cooling from free-free emission (13.09), Lr 
is the cooling from electron recombination (13.10), and Lc is the 
cooling due to collisional excitations (13.20). In the low density 
limit, each of these terms is proportional to electron density, so 
Ne cancels out. In this case, the nebular temperature will only 
depend on the energy of the ionizing photons and on the ionic 
abundances. At high densities, however, collisional de-excitations 
will decrease the efficiency of collisional cooling. Therefore, all 
things being equal, denser nebulae will be hotter. 





One very important point to consider is that collisions are, by far, 
the most important mechanism for nebular cooling. Naturally, the 
amount of collisional cooling is proportional to the metal abun- 
dance, Nz. But note: if Nz is decreased, then the temperature 
of the nebula will rise considerably. Now recall that the equation 
for collisional excitation is 





OG. 4 
qij = 8.629 x 10~° n E OR/EE gmn a (13.16) 

wile 
so if T rises, the exponential term e~4#/*7 will also rise, and 


therefore q;; will rise. As a result, it is possible that decreasing the 
number of ions will actually increase the amount of ionic emission 
from the nebula. 


Collisional cooling is always most efficient for ions where AE x 
kT. If AE > kT, then very few electrons will be excited, so 
little cooling will occur; if AE < kT, then the amount of energy 
released in the downward transition will be inconsequential. 
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The Emitted Spectrum of a Stromgren Sphere 


There are four types of emission from a Stromgren sphere 
Strong recombination lines of hydrogen and helium. 


Very strong collisionally excited lines of metals. In the optical, 
these are mostly forbidden lines of oxygen, nitrogen, sulfur, argon, 
and neon. Lines of carbon, magnesium, and silicon appear in the 


UV. 


There is some weak continuum emission associated with nebulae, 
due to free-free emission, free-bound emission, and two-photon 
emission. (The latter occurs when a 2s electron of hydrogen or 
helium decays to the 1s state by temporarily making an interme- 
diate state for itself.) At most wavelengths, this is much, much 
fainter than the line emission (although, at radio wavelengths, the 
free-free and free-bound emission is significant). 


Other weak emission lines due to resonance-fluorescence (where a 
recombination line of hydrogen happens to have the same wave- 
length of as a transition to the ground state of a metal), recom- 
bination of metals, and other odd effects. 


Note that the size of a Stromgren Sphere can be defined in two 
ways. In a radiation bounded nebula, all the ionizing photons 
are absorbed by hydrogen. Alternatively one can have a mat- 
ter bounded nebula, in which some ionizing photons escape into 
space simply because there are no hydrogen atoms around for 
them to ionize. 
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The Hydrogen Spectrum 


Since the A values for hydrogen are A ~ 10° sec~', while the q 
values are 10~*, collisional excitations for hydrogen only become 
important when 


A 
Nex — F ~ 101° em? (14.01) 


Thus, with slight exceptions for electrons in extremely high levels 
(n ~ 100), collisions do not affect the hydrogen atom. One can 
therefore write the detailed balance equation for the level popu- 
lations of hydrogen as 





CO n—1 
Np Neant (T) ae DD NO Nn t An L!,NL = Net >. N Antn'L” 


n'>n L n’=1 L 

(14.02) 
where @nz(T) is the recombination coefficient representing the 
rate at which recombinations occur directly into state nL. In 
other words, the number of recombinations directly into a level 
plus the number of radiative decays into that level is equal to the 
number of radiative decays out of the level. If (14.02) is correct, 
then the emission produced by any transition is 





=" nn! 
ea 29 NO NatAntn'L! (14.03) 
L L'=L+1 


where we have explicitly kept the individual L states separate, 
although they are degenerate in energy (and therefore summed 
for the total emission). 





Equation (14.02) represents a set of n x L linear equations with 
n x L unknowns. In theory, one can choose an electron temper- 
ature and density, look up the values of a,7 appropriate for the 
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temperature, and then solve for the population levels with lin- 
ear algebra. Once the level populations are known, the emission 
strength of each hydrogen line follow through (14.03). 





In practice, the level populations are only weakly dependent on 
Ne and Te (i.e., direct recombination is only one of N-ways into 
n = 2. Consequently, the level populations, and therefore the line 
ratios are almost entirely determined by the Einstein A values, 
not Ne or Te. 





Equation (14.02) explicitly neglects two effects. The first is col- 
lisional excitations and de-excitations. For the ground state of 
hydrogen, AE < kT, so free electrons don’t have enough energy 
to cause a collisional excitation. And, while electrons in excited 
states can be affected by collisions, the critical density for this is 
N.x 10° cm~!°. Hence, in most cases, collisions can be neglected. 
(Actually, this is not strictly true: because it doesn’t take much 
of a collision to move electrons from say, n=101 to n=100, in the 
electrons at these high levels can, and will, be redistributed. The 
result will be more like a Boltzmann distribution for these levels.) 


The other piece of physics missing from (14.02) is input to levels 
via absorptions from lower levels. Since hydrogen electrons decay 
to the ground state in ~ 107° sec™t, virtually no absorptions 
occur from these states. However, absorptions from the ground 
state can and do occur. 
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To understand the effects of absorption, let’s first consider the 
cross section of the hydrogen atom to Lyman photons. 


Line Wavelength (A) A (sect) ag (em?)  To/T 


912 A 
Ly a 1215.67 6.26 x 108 5.90 x 10714 9366 
Ly B 1025.72 1.67 x 108 9.46 x 107+5 1501 
Ly y 972.54 6.82 x 107 3.29 x 10715 522 
Ly 10 920.96 ADI x 10° 1:72 10" 27 
Ly 15 915.82 1.24 x 106 5.00 x 1071!" 8 
Ly 20 914.04 SAA x10 210x07 3 


Especially for the lower Lyman lines, the cross section for absorp- 
tion is much greater than it is for an ionization. (Recall that at 
912 Å, the ionization cross section is 6.30 x 10718 cm?.) So, if the 
on-the-spot approximation is valid, Ly absorptions must also be 
important. 


Traditionally, two types of approximations are used. Nebulae are 
divided into two types: 


CASE A: all Lyman line photons produced by recombining hydro- 
gen atoms escape the nebula. In this case, no Lyman absorptions 
occur, and equation (14.02) is strictly valid. This is the optically 
thin case. 


CASE B: all Lyman line photons are re-absorbed by other hydro- 
gen atoms. In this case, all downward transitions to the ground 
state don’t count (in the same way as ground state recombinations 
don’t count in the on-the-spot approximation). In this optically 
thick case, equation (14.02) is still valid, except transitions to the 
n=1 level are not included in the summations (since, in effect, 
they never occur). 
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In real life, almost everything is CASE B (or close to it), at least 
for the lower Lyman lines. 


CASE B has an interesting implication for astrophysics. Accord- 
ing to ionization balance, every ionization results in a recombina- 
tion (with recombinations to the ground level not counted). Once 
in the hydrogen atom, the electrons decays downward, but decays 
to the ground state don’t count. As a result, every decay must 
eventually go to the n = 2 state (since it can’t go to the n = 1 
state). All these transitions are in the optical! Thus, the net re- 
sult is every ionization produces an optical Balmer line photon. 
By counting the number of Balmer photons, you can count the 
number of ionizing photons emitted from the central star. 


The equations are even more elegant. The level populations given 
by (14.02) depend almost exclusively on the A values of the tran- 
sitions, and are only weakly dependent to Ne and Te. Since the 
level populations are (more-or-less) fixed by atomic physics, the 
line ratios are also (more-or-less) fixed by atomic physics. Thus, 
if you measure one Balmer line, you more-or-less know all the 
Balmer line strengths. So, by measuring one Balmer line (say, 
Hf), you can derive the number of ionizing photons coming from 
the central star. 
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Below is a table giving the line ratios for the lowest level Balmer 
transitions. 


Temperature 9000 10,000 20,000 

Ne (cm73) 104 102 10° 10? 104 
Ofte 5.44 3.02 3.07 1.61 1.61 
I(Ha)/I(HB8) 3.00 2.86 2.81 275 274 


HB) 0.460 0.468 0.471 0.475 0.476 
T(He)/I(H8) 0.155 0.159 0.163 0.163 0.163 


Ofte in units of cm? s~!. 


In addition to line ratios, the table also has a value called affa. 
Since the recombination coefficient for hydrogen is known, and 
since all the line ratios for hydrogen are known, it is a relatively 
straightforward task to combine these numbers and figure out how 
many recombinations eventually create an HG photon. This is 
called the “effective” recombination coefficient for HG. In general, 








AT Inn! 





n—-1 
eff > ) 
NeNpQnn/ = Notni nD = 
L=0 L'=L+1 


(14.04) 


hilar 
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Continuum Emission 





There are three sources of continuum emission from a nebula 


Free-free emission. Most of this light emerges in the radio, al- 
though some optical continuum is produced. 


Free-bound transitions (7.e., recombinations). Recall that most re- 
combinations involve slowly moving electrons, t.e., electrons with 
energies just about 13.6 eV. Recall also that there are many more 
high n levels for hydrogen than low n levels. As a result, most of 
this radition comes out at long wavelengths, though just to the 
blue of the ionization edge of each level, there is a large jump in 
the amount of continuum flux emitted. 


2-photon emission. This mechanism comes from relativistic quan- 
tum mechanics. Recall that the Hamiltonion operator for an elec- 


tron 1s 2 
Pp 


2Me 


where p is the electron momentum and ¢ is the scaler potential in 
the atom. Relativistically, however, the momentum is given by 





H = + (2, t) (14.05) 


P(x, t) (14.06) 


where ® is the vector potential of the radiation field, and q = —e 
the electron charge. Thus, the Hamiltonian operator becomes 











12 
H = |p+ zj — ed (14.07) 
or 
= (Ip? - ep) +- (P 8+8- p) + (6-8) 
~ Lame j 2MeC 2Mec2 
(14.08) 
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The first term describes the atom in equilbrium, and the sec- 
ond term governs interactions between the atom and a radiation 
field. The third, quadratic, term describes interactions between 
the atom and two photons at once. It is much smaller than the 
second term (by about 10 orders of magnitude), but this is the 
only process by which an electron in the 2s state of hydrogen can 
radiatively decay down to the 1s state. 


If the difference between the n=1 and n=2 states of hydrogen 
is 10.2 eV, then, by definition, the two photon emission process 
must create must as many photons shortward of 5.1 eV (2431 A) 
as longward. Most of this emission is therefore in the ultravio- 
let. Note, however, that since the A value for this transition is 
8.23 sec™ t}, in higher-density nebulae, electrons can be collided out 
of the 2s state before they can radiatively decay via this process. 
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Lya and Resonance Fluorescence 


Under ionization balance, every ionization is balanced by a recom- 
bination, and under Case B, every hydrogen recombination even- 
tually results in a Balmer photon (which leaves the nebula). These 
n=2 electrons will, of course, decay to n=1, but the Lya photons 
that are produced are trapped — they just get re-absorbed (and 
re-emitted) by other neutral atoms. This can happen hundreds or 
even thousands of times, but eventually, something has to happen 
to all these Lya photons. 








Actually, there are several possible ways to get rid of Lya. The 
principle ones are 


Lya may “leak” out of the nebula. These photons can random 
walk their way to the edge of the nebula, then escape, or have 
their wavelength shifted slightly due to the quantum mechanics 
of the hydrogen atom (i.e., they can be emitted “in the wings” of 
the line). The process is assisted if the nebula has bulk motions, 
which can redshift a Lya photon out of the line-center. 





Lya may get destroyed by hitting a dust grain. There’s usually 
some dust floating around, and if a photon hits it, the photon will 
be destroyed (and the grain will be heated). One dust grain can 
kill many, many Lya photons. 


The electron may be collided out of the n=2 P-state into either 
the n=1 state, or the n=2 S-state. In the latter case, the electron 
may (if it doesn’t get hit again) decay via 2-photon emission. 


A “resonance fluorescence” occurs. Actually, this doesn’t happen 
with Lya per se, but with Lyg (destroying the photon before it 
gets to n=2) or with the Lya of He II. (Helium has the same 
problem as hydrogen, and has its own Case B set of equations.) 
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The way resonance fluorescence works is this. Let’s take He II Lya 
as an example. The wavelength of this Lya transition is 303.78 A. 
A permitted, ground-state transition of OT* (3d ° P) —> 2p? 3 P2) 
is 303.80 A. This coincidence can cause the helium Lya photon to 
excite oxygen instead of a helium ion. The oxygen atom can then 
decay (via another route) and emit photons of various wavelengths 
that escape the nebula. (There are a whole bunch of permitted 
oxygen emission lines between 2800 A and 3800 A that are pro- 
duced in this way.) 


Resonance fluorescence is more interesting than it is useful. Unless 


one is doing detailed radiative transfer calculations, it can usually 
be ignored. 
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Nebular Diagnostics 


Ions with ground state configurations of p?, p3, and pt have four 
low-lying levels with energy separations of the order of kT. The 
p? and pt configurations are especially good for determining elec- 
tron temperatures; ions with this configuration include Nt, OTT, 
Net4, St+, and Art‘ (all p?) and O, Nett, and Art (p*). Ions 
with p? configurations, such as Ot, Nett+, St, and Artt+, are 
useful for measuring density. 


The flux generated by any forbidden transition (from i to 7) is 


where N; is the number of ions in the 7 state, A;; is the Einstein 
A value of the transition, and hvi; the energy difference of the 
transition. Also recall that the number of ions in state 7 is given 
by the conditions of detailed balance 


Ss” NeNjQ5i + x. NV Ag = ` NeNiqij + ye NA (13.28) 
j#Ft j>i j#Ft j<i 
In other words, the number of collisions into the level plus the 
number of radiative decays into the level is equal to the number 
of collisions out of the level plus the number of decays out of the 
level. The collisional rates, qj; are given by 
Oli; 
(3) cm? s" fori>j (13.15) 
1/2 
wiTe 





qij = 8.629 x 107° 





and 





E 
qij = 8.629 x 10~° OH) e-AB/kT cms! fori <j (13.16) 
Wile 


Thus, for any Te and Ne, one can solve for N;, and thus for the 
F;j. By taking the ratios of certain lines, measurements of T} and 
Ne can be made without knowing distance or abundance. 
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Estimating the Electron Temperature 


Shells with p? or pt configurations have a 3P term for the ground 
state, a singlet D term AE ~ kT above the ground state, and a 
singlet S term with AE ~ kT above the D state. This type of 
3-1-1 splitting works well for temperature determinations. 


Consider the OTT ion. The A5007 
line and the 44959 both start at the 
1 Də level, so the ratio of these two 
lines is simply the ratio of their A 
values. Since the A value for A5007 
is 3 times that of the 4959 transi- 
tion, 1(5007)/I(4959) = 3. Similarly, 
the fraction of tSo ions that produce 
optical (A4363) light rather than UV 
(A2321) light is 


A4363 


= ——— 14.10 
A4363 + A2321 ( ) 





A 4363 








A 5007 





À 4959 
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p< 





(Actually, the denominator of equation (15.10) should also include 
the A values to the 3P) and *P states, but these rates are negli- 


gible compared to the other two.) 





Now the ratio of ?P —! D collisions to °P —-! S collisions is 


R= 


T1/2 


_ OCPD) AEip15/kT 
QP," S) 
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Noar (ERPE) BE Enas 


(14.11) 


If the number of radiative decays into 'Dz from tSo is small com- 
pared to the number of collisions into t Də from the ground state, 
then the intensity of each line is simply proportional to the num- 
ber of collisions into the level. So if we combine (14.10) with 
(14.11), then 





T(\4959) + 1(A5007) Q(@P,! D) A4959 + Asoo7 7 AE/kT 
T(A4363) OG PLS) Aiea Vises. 
(14.12) 
where AF is the energy difference between the | Dz state and the 
'So state, and Ð is the “mean” frequency of a photon produced 


by a decay from ! Da, i.e., 





A A 
i V4959A4959 + 1500745007 (14.13) 
A4959 + As5007 


This expression, of course, assumes that all the D and S electrons 
come from collisions with ground state electrons. It therefore is 
strictly true only in the low-density limit. At higher densities, 


(Ne 10°), the line ratios must be solved via the detailed balance 
equations. 
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Temperature T 
The [O III] lines intensity ratio (in the low density limit) as a 
function of temperature. Note that at low temperatures, very 
few electrons are collisionally excited up to S-term level. Conse- 
quently, the 44363 line is extremely weak. Similar relations exist 
for other ions with p? or p* configurations. 
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Estimating the Electron Density 


Ions with p? configurations have a pair of closely spaced 7D levels 
AE ~ kT above the ground state, and another pair of 7P levels 
AE ~ kT above the D term. We can use the 1-2-2 splitting to 
measure the electron density. 








To see how, consider the OT ion. The 12 2P 
2D term is split into two levels with 

essentially the same energy above the 

ground state. The ° D3 j2 level has a 

statistical weight of w = (2x $)+1 = 4; 

the *Ds/2 level has w = 6. Since the 52 "D 
ground state of a p3 configuration is 
an S state (i.e., L = 0), the collision 
strength from the ground state is the 13729 13726 
same for both levels, except for the sta- 
tistical weight (see equation 13.18). 
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So, if Q is the total collision strength from S to D, the collision 
strength from S' to *Ds/2 is 0.69, and that from S' to *Ds/2 is 
0.49. (In other words, according to the definition of statistical 
weights, six-tenths of the D states have J = 5/2, and four-tenths 
have J = 3/2: the electrons don’t care which state they enter, 
and there are more J = 5/2 states than J = 3/2 states.) 





Now, in the low density limit, every collision upward results in a 
radiative decay downward. The ratio of ? D; /2 9? S3 /2 decays to 
*Ds /2 a5 Ss 2 is then simply the ratio of the number of collisions 
into each level. Since the energies of both states are essentially 
identical, the ratio of the lines is just the ratios of Q, which are, 
in turn, just the ratios of the statistical weights. Therefore, in the 
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low density limit, 


I(A3729) W3729 6 
a Se ee 14.14 
I(A3726) W3726 4 ( ) 





In the high density limit, however, things are different: elec- 
tron collisions are continually populating and depopulating the 
D states. Again, since the statistical weight of ? D; J2 is 1.5 times 
that of 7 Ds /2, all things being equal, there will be 1.5 more decays 
from that level. However, since each level has an infinite supply of 
electrons from collisions, the ratio of the lines also depends on the 
ratio of the A values (see equations 14.09 and 14.10). (Imagine if 
one level had an A value ten times that of the other; in that case, 
10 photons would be produced in the time it takes the other level 
to produce 1 photon.) Thus, in the high density limit 





I(A3729) = W3729A3729 =| p 43729 
I (A3726) W3726A3726 l A3726 





(14.15) 





In the intermediate case, of course, one must solve for the line 
ratios via the equations of detailed balance. However, the ratio 
of the two lines is smoothly function, so simple interpolations will 
work. 
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Electron density (cm~?) 
The variation of the intensity ratio of [O II] 43729 to 43726, and 
of [S II] 46716 to 6731 as a function of density at T = 10,000 K. 
At other temperatures, the curves are very nearly identical, if the 


x-axis is defined as Ne/4/T/10, 000. 
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Measuring the Ionizing Flux of Photons 


In the optical thick CASE B is valid, then it is easy to use the 
recombination lines to measure the ionizing flux from the central 
star (or stars). Recall that from ionization equilibrium 


oO ies R 
/ z = Q(H’) = f Ne Npa gdV (12.36) 
V 0 


(0) 


where V is the nebular volume. (In the spherical nebula case, 
dV = 4rr?dr.) Also recall that under CASE B, every ioniza- 
tion causes a recombination which creates a Balmer photon that 
escapes the nebula. The H8 luminosity is therefore 


R R 
L(HB) = / AnjygdV = hvyg i NeNpoyig dV (14.16) 
The ratio of ionizing photons to H8 photons is therefore 


R e 
LHS) /hvas _ Jo NeNpoyis aV 


14.17) 
L R ( 
Loret fy NeNpapdV 
In other words, the ratio of HG photons to ionizing photons is the 
density weighted average of ae divided by the density weighted 
average of ag. So 


L(HB)/hivng _ og 


oo L, 
toe not AB 


(14.18) 


Both an and ap are physical functions that depend only on 
temperature (in a similar fashion). Thus, given an estimate for 
the nebular temperature, the constant relating HG photons to 
ionizing photons is easily found. 


51 


Measuring the Temperature of an Ionizing Source 


If CASE B holds, and the ionizing source is a single star, we can 
use the nebulas H£ flux and the V-band magnitude of the cen- 
tral star to measure the central stars temperature. First, let’s 
approximate the flux distribution of the ionizing star as a black- 
body, B(T). The ratio of V-band photons to ionizing photons 
will clearly be a function of temperature. (You can compute this 
on a calculator, if you like.) Let’s call that function F(T). From 
(14.18) 


F(T) = A a (14.19) 


eBay (T saan o> aid 
/ v ( m I ae Y Pu “He 
w M w IW A Cree 


Ly ap 


PD = CURB) asf 














or 


hvus (14.20) 


We can easily convert total luminosity to observed flux by dividing 
both numerator and denominator by distance squared. Therefore 


fv aB 
POL) = FHA) a cil hig (14.21) 
In theory, the observed V-band flux of the central star and the 
observed H 8 flux from the nebula are both visible. One can there- 
fore find the value of the right-hand side of the equation, and look 
up the temperature of the source via the function F(T). This is 
called the Zanstra method for measuring stellar temperatures. 


Variants of the Zanstra method work for He I and He II as well. 


(He I has a line at A5876 which serves the same function as HØ; 
the He II equivalent of H£ is a line at A4686.) 


52 


Estimating Extinction to a Nebula 


Estimating the extinction to a nebula is almost trivial. The ra- 
tio of the Balmer emission lines is extremely insensitive to den- 
sity and temperature. You cannot go to far wrong if you assume 
Ha/H6 ~ 2.86. If the ratio is larger than this, the discrepancy 
is almost certainly due to differential extinction, t.e., reddening. 
Recall that the amount of extinction (in magnitudes or log flux) 
at any (optical) wavelength is A, œ 1/X. According to this 1/A 
extinction law P 

= log = A (14.22) 


obs 


log fX 


where qo represents the observed flux at any wavelength, fy 
is the de-reddened flux, and a is a constant that represents the 
amount of extinction. Now, for consistency with the literature, let 
c equal the total logarithmic extinction at H8, i.e., c = a/Aage1.- 


So 





log fab = log fag —C (14.23) 


In terms of c, the ratio of the observed flux at any wavelength A 
to the observed flux of HG is then 


log it = log abs =c (> — 1) (14.24) 
F fag AHB 


Since the intrinsic Ha/Hf ratio is known, (14.24) can immediately 
be solved for c. Once c is known, the de-reddened H flux fol- 
lows immediately from (14.23), and, with c and fy g both known, 
(14.24) can be used to solve for the de-reddened flux at any wave- 
length. 
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Deriving Nebular Abundances 


The key to deriving nebular abundances is to use the nebular 
diagnostics to fix the electron temperatures and densities. Once 
that is done, the abundance of any species follows simply from the 
observed fluxes. For instance, for hydrogen, He I, and He II, the 
amount of emission is 


I(H8) = NeNp aff > hvug (14.25) 
1(\5876) = NeNge+ aff- - hvsere (14.26) 
T(A4686) = NeNize++ OShbeg < hvasse (14.27) 


while for collisionally excited ions, the amount of emission is 
dj; = NiNeqij $ hvi,j (14.28) 


Since the values of a and q;; are functions of temperature and den- 
sity only, once those values are known, the line strengths become 
directly proportional to abundance. 


Note that if one wants to measure the total abundance of oxygen, 
one needs to add up the abundances of oxygen in all its states, 
i.e., O, OT, OTF, etc. Fortunately, by observing the strengths of 
the recombinational lines of hydrogen, He I, and He II one can 
get a pretty good idea of the underlying distribution of ionizing 
photons. Thus, one can estimate the relative abundance of, say, 
Ot? to Ot?, even if no lines are OT? are observable. 
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One problem that arises in nebular analysis is that there is quite 
a lot of redundancy in the information provided by the emission 
lines. In some cases, not all the measurements of density and 
temperature are consist. In particular, the density measurements 
from the [O II] lines are sometimes much larger than those mea- 
sured via other techniques (specifically, radio observations of the 
nebular continuum). This has lead to the introduction of f, the 
filling factor. Many nebulae are therefore modeled by assuming 
the existence of clumps of material with density Ne separated by 
regions with density zero. This value f finds its way into many of 
the equations presented, i.e., the equation of ionization balance 
becomes 








Q(H°) = ST RSNN po -f (14.29) 
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Summary of Emission Line Physics 


e The recombination lines of hydrogen and helium allow you to 
count the number of photons shortward of their ionization edge 
and therefore map out an objects flux distribution in the far-UV. 
The interpretation of these lines, does, however, depend on the 
density squared. 





e The Balmer (and helium) line ratios are (to a very high degree) 
fixed by atomic physics, and do not depend on nebular conditions. 
In particular Ha/H@ = 2.86 (with at most a couple percent er- 
ror). If one measures a higher ratio, it is almost certainly due to 
interstellar reddening. 





e The collisionally excited forbidden lines of oxygen, nitrogen, 
sulfur, argon, and neon allow one to measure nebular conditions. 
The 2-2-1 structure of the p? configuration is particular sensitive 
to the density of free electrons. The ratio of [O II] 43729/A3726 is 
best for this, although because the lines are so close together, they 
are often blended into one ‘3727’ line. [S II] 6717/6731 is just 
as sensitive, but may be difficult to measure since the abundance 
of sulfur is generally just ~ 4% that of oxygen. 








Ions with p? and p* configurations have a 3-1-1 structure, and 
therefore can be used to measure the electron temperature of the 
ISM. Typically, this is done with OTT via the ratio of A5007 + 
A4959 to A4363. Note, however, that the 44363 line will be be- 
tween 30 and 1200 times weaker than 5007. So its detection can 
be challenging. (NT can also be used with 45755 as the key line, 
but often times, most of the nitrogen is doubly, rather than singly 
ionized. ) 


e If one can fix the electron temperature and (to a lesser ex- 
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tent) the density, then the ionic abundaces follow directly from 
the strength of the forbidden line. There is still the issue of how 
much gas is in the various ionization states, but if one measures 
abundances for two species (say, OT and OTF) one can usually 
solve this problem. 


e If one cannot fix the electron temperature and density, then 
a direct measurement of abundance is impossible. Nevertheless, 
there are a number of empirical relations between abundance and 
strong-line strengths, such as the Rə3 relation 


I({O 11]A3727) + I({O 111]A5007) 


a THA) 





These have been derived locally using various samples of ionized 
regions. The applicability of these relations to specific problems 
may be controversial. 
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e The ratios of strong lines have also (empirically) been shown to 
be useful for object classification. This was first done by Balwin, 
Phillips, & Terlevich (1981); H II regions, planetary nebulae, ac- 
tive galactic nuclei, and LINERS fall on different parts of these 
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“BPT” diagrams. Note that the line ratios used are quite close 
together in wavelength; thus, each index is independent of red- 
dening. Again, these relations are completely empirical, and thus 
their applicability to certain objects may be challenged. 
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